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Sum ma ry 

A conical Euler code was developed to study unsteady vortex-dominated flows 
about rolling highly-swept delta wings, undergoing either forced or free-to-roll 
motions including active roll suppression. The flow solver of the code involves a 
multistage Runge-Kutta time-stepping scheme which uses a finite-volume spatial 
discretization of the Euler equations on an unstructured grid of triangles. The code 
allows for the additional analysis of the free-to-roll case, by including the rigid-body 
equation of motion for its simultaneous time integration with the governing flow 
equations. Results are presented for a 75° swept sharp-leading-edge delta wing at a 
freestream Mach number of 1.2 and at a = 10° and 30° angle of attack. At the lower 
angle of attack of a = 10°, a forced harmonic analysis indicates that the rolling moment 
coefficient provides a positive damping which is verified in a free-to-roll calculation. 
In contrast, at the higher angle of attack of a = 30°, a forced harmonic analysis indicates 
that the rolling moment coefficient provides a negative damping at the small roll 
amplitudes. A free-to-roll calculation for this case produces an initially divergent 
response, but as the amplitude of motion grows with time, the response transitions to a 
wing-rock type of limit cycle oscillation. The wing rocking motion may be actively 
suppressed, however, through the use of a rate-feedback control law and 
antisymmetrically-deflected leading-edge flaps. The paper provides descriptions of the 
conical Euler flow solver and the free-to-roll analysis. Results are presented which 
give insight into the flow physics associated with unsteady vortical flows about forced 
and free-to-roll delta wings, including the active roll suppression of this wing-rock 
phenomenon. 



Introduction 


In recent years, the understanding and prediction of the complex flows about 
modern aircraft at high angles of attack have been research topics that have generated 
much interest within the fluid dynamics community. 1 >2 These aircraft typically have 
thin highly-swept lifting surfaces which produce a vortical flow over the leeward-side 
of the vehicle at high angles of attack. This vortical flow can have beneficial effects on 
performance, such as lift augmentation at high-a, but may also have adverse effects 
such as structural fatigue due to tail buffet and also stability and control problems such 
as wing rock, wing drop, nose slice, and pitch-up. 8 Consequently, considerable work 
has been done experimentally to try to understand the basic flow physics involved in 
vortical flows about delta wings at high angles of attack. Experimental research efforts 
directed towards understanding and documenting steady vortical flows are typified by the 
detailed flowfield measurements about simple, cranked, and canard-wing delta 
configurations at low speed by Hummel 4 and low speed tests on a 75° swept delta wing by 
Kjelgaard and Sellers. 5 For supersonic-freestream Mach numbers, vortical flows have 
been measured by Squire 5 for an elliptic cone delta wing and by Miller and Wood 7 for a 
series of swept sharp-leading-edge delta wings. Efforts on investigating unsteady 
vortical flows experimentally have been reported by Nguyen et al. 8 for forced harmonic 
and free-to-roll motions of an 80° swept delta wing in low speed flow. In Ref. 8, the 
wing was found to undergo self-induced periodic roll oscillations known as wing rock, 
for angles of attack greater than 25°. Levin and Katz 8 tested both 76° and 80° swept 
delta wings and found that only the 80° model would exhibit wing rock at high-a. 
Further studies have been performed by Nelson and co-workers 10 - 11 at Notre Dame 
University. These studies have shown, for example, the time histories of the vortex core 
position during a cycle of wing rock 1 0 and the static and dynamic effects due to vortex 
breakdown. 11 Also Ng et al. 18 have recently reported experimental results obtained in 
a water tunnel which show wing rock for several different delta wing planforms along 
with detailed flow visualization diagrams. These studies 4 ’ 18 have contributed 
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significantly to the understanding of steady and unsteady vortex-dominated flowfields, 
although much work remains to be done. 

From a computational point of view, considerable effort has also been spent on 
developing methods of predicting steady and unsteady vortex-dominated flows. 1 ®- 14 
Hoeijmakers, 1 3 for example, gives a review of computational methods for the 
determination of steady vortical flow characteristics with an emphasis on classical 
methods including discrete vortex, cloud-in-cell, panel, vortex layer with finite core, 
leading-edge suction analogy, and vortex-lattice. With respect to unsteady methods, Hsu 
and Lan 15 presented a nonlinear mathematical model for calculating wing-rock 
characteristics based on aerodynamic derivatives evaluated using steady-flow 
aerodynamics at average dynamic conditions. Researchers at Virginia Polytechnic 
Institute and State University 1 ®' 1 ® have simulated wing rock using an unsteady 
vortex-lattice method to predict the aerodynamic loads and have integrated the equation 
of rolling motion using a predictor-corrector method. The methods of both Ref. 15 and 
Refs. 16-18 predicted with reasonable accuracy the low-speed wing rock 
characteristics of the delta wings studied in Refs. 8 and 9. Use of the more modern 
computational fluid dynamics (CFD) techniques for the prediction of vortex-dominated 
flows 14 have primarily focused on steady applications, 1 9 ' 2 ® although there are 
notable exceptions where applications have been made to rolling delta wings undergoing 
forced harmonic 27 ' 29 and free-to-roll® 9 motions. Kandil and Chuang, for example, 
have reported results for rolling delta wings obtained using the conical Euler equations 
for sharp-leading-edge wings 27 and the conical Navier-Stokes equations for rounded- 
leading-edge wings. 2 ® Batina 29 has also presented results for a rolling delta wing 
obtained using a conical Euler flow solver based on the use of unstructured grids of 
triangles. Lee and Batina® 9 extended the methods of Ref. 29 to include a free-to-roll 
capability and showed results for a freely rolling delta wing that exhibited a limit cycle 
or wing-rock type motion. 
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The objective of the current research is to study unsteady vortex-dominated 
flowfields by using the conical Euler equations as an efficient first step to investigating 
the full three-dimensional problem. Reference 30 reported the development of a conical 
Euler code to study unsteady vortex-dominated flows about rolling highly-swept delta 
wings undergoing forced harmonic and free-to-roll motions. The purpose of this paper 
is to report the extension of these methods for the study of flows about delta wings 
undergoing pulsed motions and free-to-roll motions including roll suppression. The flow 
solver of the code is similar to that of Ref. 29, which involves a multistage Runge-Kutta 
time-stepping scheme and a finite-volume spatial discretization of the Euler equations 
on an unstructured grid of triangles. The code was modified to allow for the additional 
analysis of the free-to-roll case, by including the rigid body equation of motion for its 
simultaneous time integration with the governing flow equations. The analysis also 
includes a capability for implementing an active feedback control law with 
antisymmetrically-deflected leading-edge flaps for roll suppression of the wing- 
rocking motion. Although limited experimental work has been conducted on the use of 
flaps 31 and leading edge blowing 32 - 33 for roll control, even less analytical and 
numerical work 34 has been performed. Results are presented here for a 75° swept 
sharp-leading-edge delta wing at a freestream Mach number of 1.2 and at a = 10° and 
30° angle of attack. The results demonstrate the successful CFD simulation of a wing- 
rock phenomenon including its active roll suppression. The paper gives descriptions of 
the conical Euler flow solver and free-to-roll analysis. The results give insight into the 
flow physics associated with unsteady vortical flows about forced and free-to-roll delta 
wings, including the active roll suppression of this wing-rock phenomenon. 

Symbols 

a„ freestream speed of sound 
rolling moment coefficient 

C^ g rolling moment coefficient transfer function due to flap deflection 
4 



rolling moment coefficient transfer function due to roll 

C root chord of wing 

e total energy 

i imaginary unit, satisfying the equation i 2 = -1 

I xx mass moment of inertia about longitudinal axis 


k 

K v 

l 

M k 

q~ 

s 


reduced frequency based on one half of the root chord, 


coc 


2U 


oo 


control gain 


rolling moment, positive clockwise when viewed from aft 
freestream Mach number 
freestream dynamic pressure 

planform area 


t nondimensional time, — — 

c 

U component of velocity vector in x-direction 

v component of velocity vector in y-direction 

w component of velocity vector in z-direction 

a angle of attack 

Y ratio of specific heats 

8 leading-edge flap deflection angle, positive clockwise when viewed from 
aft 

At nondimensional time step 

p. x structural damping 

p density 

freestream density 

(J) instantaneous roll angle, positive clockwise when viewed from aft 
<|) 0 harmonic and pulse roll angle amplitude 
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Governing Equations 


The flow is governed by the time-dependent Euler equations which may be written in 
conservation law form as 


3Q 3E d F 3G n 
+ =0 

where Q is the vector of conserved variables defined by 

Q = [p ( pu,pv 1 pw,e] T 


(1 ) 


( 2 ) 


and E, F, and G are the convective or inviscid fluxes given by 


pU 

pUu -i-p 
E=l pUv 
pUw 

(e +p)U + x t p 

pV 
pVu 

F=| pVv+p 
pVw 

(e+p)V+y t p 
pW 

pWu 

G= pWv 
pWw +p 
[(e +p)W +z t pj 


The contravariant velocities U, V, and W are defined by 


U = u — x t V = v - y t W = w-Z{ 


(3a) 


(3b) 


(3c) 


(4) 
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where x t ,y t and z t are the grid speeds in the x, y, and z directions, respectively. The 
pressure p is determined by the equation of state for a perfect gas 


p = (y- 1) [e - j p (u 2 + v 2 + w 2 )] 


(5) 


and the equations have been nondimensionalized by the freestream density and the 
freestream speed of sound. 

If interest is restricted to supersonic flow past conical bodies, then the conical flow 
assumption can be made. This reduces the problem from three dimensions to two 
dimensions, which significantly decreases the computational resources that are required 
to investigate such flows. The conical flow assumption is exact for steady inviscid 
supersonic flow. For unsteady flows, however, the conical assumption implies an 
instantaneous propagation of disturbances in radial directions. The conical flow 
assumption involves a change of variables according to 


t y „ z 

^ = x 11 = - C=- 

X X 

The three-dimensional Euler equations then reduce to 

ff*3T (F -’ lE) + ic (G - CE)+2E '° 


( 6 ) 


(7) 


Equation (7) may be rewritten in integral form for solution as 
y f Q dll dC + J [( F - r,E) dC - ( G - C E) dll ] + J 2 Edt, dC = o ( 8 > 

‘a a 

where the second integral is a boundary integral resulting from application of the 
divergence theorem. 
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Euler Solution Algorithm 


The unsteady conical Euler equations are solved using the multistage Runge-Kutta 
time-stepping scheme of Ref. 29. This algorithm uses a finite-volume spatial 
discretization for solution on an unstructured grid made up of triangles. The original 
algorithm of Ref. 29 was a node-based scheme whereby the flow variables are stored at 
the vertices of the triangles. A second algorithm, a cell-centered scheme, was employed 
in the present study. This second scheme is based on unpublished work of the second 
author. In the cell-centered scheme, the flow variables are stored at the centroids of the 
triangles. In both algorithms, artificial dissipation is added explicitly to prevent 
oscillations near shock waves and to damp high-frequency uncoupled error modes. 
Specifically, an adaptive blend of harmonic and biharmonic operators is used, 
corresponding to second and fourth difference dissipation, respectively. The biharmonic 
operator provides a background dissipation to damp high frequency errors and the 
harmonic operator prevents oscillations near shock waves. The algorithms also employ 
enthalpy damping, local time stepping, and implicit residual smoothing to accelerate 
convergence to steady state. The local time stepping uses the maximum allowable step 
size at each grid point for the node-based scheme and for each triangle in the cell- 
centered scheme, as determined by a local stability analysis. The implicit residual 
smoothing permits the use of local time steps that are larger than those imposed by the 
Courant-Friedrichs-Lewy stability condition. This is achieved by averaging the 
residuals implicitly with neighboring values. A time-accurate version of the residual 
smoothing is also used for global time-stepping during unsteady applications of the code. 

With respect to boundary conditions, freestream conditions are applied along the 
farfield boundary, and a reasonably large computational grid is used so that the bow 
shock is captured as a part of the solution. A flow tangency (or slip) condition is applied 
to the inner boundary which represents the wing. Also, for unsteady applications, the 
grid is moved to conform to the instantaneous position of the wing using a combination of 
rigid body rotation and deforming mesh movement. A deforming mesh algorithm 
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discussed below is used to deform the grid about the deflected leading-edge flaps in a 
wing-fixed coordinate system. This deformed mesh is then rotated as a rigid body to 
account for the instantaneous roll angle of the wing. In this application, grid speeds are 
computed at the nodes and are included in the governing equations to account for the 
relative motion between the grid and the fluid. 

Deforming Mesh Algorithm 

The deforming mesh algorithm as developed in Ref. 35, models the triangulated 
mesh as a spring network where each edge of a triangle represents a spring with 
stiffness inversely proportional to the square of its length. In this method, the grid 
points along the outer boundary are held fixed while the grid points along the wing 
(inner boundary) are specified. The locations of the interior points are then determined 
by solving the static equilibrium equations which result from a summation of forces at 
each node in both the ti and £ directions. The solution of the equilibrium equations is 
carried out using a predictor-corrector method which first predicts the new locations of 
the interior points by extrapolation from the previous time levels and then corrects 
these locations by using several Jacobi iterations of the static equilibrium equations. 
The predictor-corrector procedure is relatively efficient since it requires only a few 
Jacobi iterations to efficiently move the mesh. 

Pulse Transfer-Function Analysis 

Generally, unsteady load coefficients can be obtained by calculating several cycles of 
a harmonically forced oscillation with the determination of the load based on the last 
cycle of oscillation. This method of harmonic oscillation requires one flowfield 
calculation for each value of reduced frequency of interest. By contrast, the unsteady 
load coefficients may be determined for a wide range of reduced frequency in a single 
flowfield calculation using the pulse transfer-function analysis. The pulse transfer- 
function analysis has been employed in the past to determine the generalized 
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aerodynamic- forces (GAFs) which are used in aeroelastic analyses. 36 - 37 in the pulse 

analysis, the unsteady force coefficient is computed indirectly from the response of the 

flowfield due to a smoothly varying, exponentially-shaped pulse. Results computed 

using the pulse analysis for a pitching flat-plate airfoil were presented in Ref. 36. 

These results were in good agreement with parallel linear theory calculations which 

validated the accuracy of the analysis. Applications to transonic airfoil cases were also 

in good agreement with the GAFs computed using the harmonic method which tends to 

verify that the analysis is valid for predicting the small perturbation response about a 

nonlinear flowfield. 37 Therefore, because of the computational efficiency of the pulse 

transfer-function analysis, the capability was implemented within the conical Euler 

code to calculate the rolling moment coefficient due to roll, C/ , of a delta wing. The 

*♦ 

pulse is expressed as 

i/T\ , -Mi(t-t c ) 2 

<xt) = 4>0 e (9) 

where <}> 0 is the pulse amplitude, is the freestream Mach number which determines 
the width of the pulse, and t c is the nondimensional time at the center of the pulse. A 
small rolling pulse is prescribed for the delta wing, and the aerodynamic transient is 
calculated. The aerodynamic transient is then used to obtain the rolling moment 
coefficient in the frequency domain by a transfer-function analysis. In this case, a fast 
Fourier transform (FFT) of the rolling moment coefficient is divided by a FFT of the 
pulsed rolling motion. The transform assumes that the system is locally linear which is 
shown to be a valid assumption for the pulse amplitude of 1° used in the present study. 

Free-to-Roll Analysis 

In this section, the roll equation of motion, the time-marching solution procedure 
and the active roll suppression are described. 

Roll Equation of Motion 


The equation of motion for a rolling delta wing can be expressed as 



Ixx 4 , — Hx 


( 10 ) 

where <t» is the roll angle which is positive clockwise when viewed from aft, l xx is the 

mass moment of inertia about the longitudinal axis, l is the aerodynamic rolling moment 
also positive clockwise, and pi x is a structural damping term (dot superscripts indicate 

differentiation with respect to time). In order to nondimensionalize Eq. (10), the 
angular rates are multiplied by the root chord of the delta wing, c, and divided by the 
freestream speed of sound, a„ . The rolling moment coefficient is defined as 


q«, Sc (11) 

where is the freestream dynamic pressure and S is the planform area. The 
nondimensional rolling equation of motion can then be written as 




-Ci C, 


C 2 4» / 


( 12 ) 


where 


Ci = 


Mi sc 3 P , 


21 


xx 


(13a) 


C 


2 “ 


HxC 

lxx 


(13b) 


Note that the prime superscripts indicate differentiation with respect to nondimensional 
time, t. The structural damping term is added to simulate a sting bearing mount. This 
type of bearing mount was used in the low-speed wind tunnel investigations of wing rock 
reported in Refs. 8-11. 


Time-Marchina Solution 

The solution procedure for the time integration of Eq. (12) is based on a finite 
difference representation of the time derivatives. The time derivatives are expressed in 
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terms of second-order-accurate finite-difference approximations. After substituting 
these expressions into Eq. (12), the roll angle at time level n+1 can be expressed in 
terms of the roll angle at previous time levels as 

<t> n+1 =[C 1 Cj + 1 At 2 + (5 + 2 C 2 At) <J> n 

-(4 + |C 2 At)*})"' 1 + <}>"' 2 ]/[§C 2 At + 2] (14) 

The rolling moment at time level n+1, C" +1 , is estimated from a linear extrapolation of 
at the previous two time levels. This predicted value of is used to determine the 
roll angle at time level n+1, <J) n+1 . The flowfield is then calculated about the wing at 
this roll angle, and the actual value of the rolling moment coefficient is determined. The 
rolling moment coefficient is then updated for use in the next time step. Due to the 
explicit time-marching of the Euler code used in this study, the time steps required for 
stability were small, and thus, it was not necessary to iterate between the roll angle 
calculation and the flowfield calculation at each time step. Previous studies of time- 
marching aeroelastic analyses using a similar explicit scheme have shown this to be the 
case (R. D. Rausch: Personal Communication, October 31, 1989). For a free-to-roll 
calculation, steady-state initial conditions are specified for <() — <J)° , CJ 1 , andC^ 0 . An 

initial angular velocity is imposed to provide an initial perturbation to the wing. 

Active Roll Suppression 

Active roll suppression is achieved through the addition of an active rate-feedback 
control law to the time-marching solution procedure. A simple control law was chosen 
of the form 

5 = K V <J>' (15) 

where K v is the control gain and 5 is both the left and right leading-edge flap deflection 
angles measured positive clockwise from the flap hinge lines. The control law is applied 
to both left and right flaps simultaneously which results in an antisymmetric 
configuration. The time-marching solution procedure is the same as that which was 


described in the preceding section. However, after the roll angle at time level n+1 is 
determined from Eq. (14), the flap deflection angle is determined from Eq. (15) using a 
second-order-accurate finite-difference expression for the angular velocity <)>'. The 
deforming mesh algorithm is then applied in addition to the rigid rotation to move the 
mesh to its new position. As before, the flowfield is calculated about the wing at its new 
position, and the rolling moment coefficient is determined and then updated for use in the 
next time step. The same initial conditions as described in the preceding section are 
applied to begin the calculation. 


Results and Discussion 

Calculations were performed for a 75° delta wing at a freestream Mach number of 
1.2 and at a = 10° and 30° angle of attack. The wing has thickness and sharp leading 
edges as indicated in the partial view of the grid shown in Fig.1. The thickness-to-span 
ratio at this cross section is 0.025 and the lower edge bevel angle is 10°. The grid, 
which was generated using an advancing front method, 38 has a total of 4226 nodes and 
8299 elements. The grid was designed to be fine on the leeward side of the wing where 
the dominant flow features are expected to occur and to be coarse on the windward side of 
the wing where the flow gradients are expected to be small. As discussed previously, for 
unsteady applications, the mesh is rotated as a rigid body to conform to the instantaneous 
position of the main part of the wing. The mesh is deformed locally near the leading 
edges to conform to the instantaneous position of the flaps. As examples of mesh 
movement, partial views of the left leading-edge flap at a positive (8= 10°) and a 
negative (8 = -10°) flap deflection angle are shown in Figs. 2(a) and 2(b), 
respectively, with the wing rotated through 10° of motion. The hinge point of the flap 
coincides with the inboard bevel edge and so the flap length is approximately 28% of the 
semi-span. As shown in the figure, the mesh moves smoothly as the wing rolls and the 
flaps are deflected. 



Steady and unsteady results including the pulse, forced harmonic and free-to-roll 
calculations are presented for both a = 10° and 30°. The rate-feedback control law is 
applied to the a = 30° case since it was the only free-to-roll case to exhibit a wing- 
rock behavior. 


Steadv-State Results 

Steady-state results were obtained to determine the basic character of the vortical 
flows and to provide starting solutions for the unsteady cases. A comparison of total 
pressure loss contours from these solutions, shown in Fig. 3, illustrates the effects of 
angle of attack. For the a = 10° angle of attack shown in Fig. 3(a), the contours indicate 
that the flow separates from each of the leading edges of the wing producing two small, 
widely-spaced circular vortices. In contrast, at the a = 30° angle of attack, shown in 
Fig. 3(b), the contours indicate that the flow separating from the leading edges produces 
two large, more closely-spaced circular vortices. Also, as the flow accelerates beneath 
the vortices of Fig. 3(b), vertically-oriented crossflow shock waves are formed on the 
outboard portions of the wing. Weaker shock waves are formed on the top of each vortex. 
These vertically-oriented shocks are located above the core of the vortex. A weaker 
horizontal shock wave is also present between the vortices. 

Pulse Transfer-Function Results 

The pulse transfer-function analysis was performed to determine the small 
amplitude stability and response characteristics of the wing. A stability analysis was 
derived by first recalling the nondimensional rolling equation of motion given by 

<J) = C 1 - C 2 (12) 

For simplicity, the structural damping term (C 2 ) is set equal to zero resulting in 



Assuming that the rolling moment coefficient can be written as the product of the rolling 

moment coefficient transfer function C / and the roll angle <(>, 

4 <t> 

C t =C t^ (17) 

then for simple harmonic motion 

C g =[Re(C^)+ilm(C^ )]<}> (18) 

where Re(C^) and Im(C^) represent the real and imaginary parts of the first 

harmonic component of C^, respectively. In this case, the real part of the rolling 

moment coefficient transfer-function represents an aerodynamic stiffness and the 
imaginary part represents an aerodynamic damping. Therefore, for this simple one 
degree-of-freedom case, the sign of Im(C^) determines the stability in roll of the 

wing for small perturbations. In other words, a negative Im(C^) indicates a positive 

aerodynamic damping which would cause a free-to-roll wing to be stable, and a positive 
Im(C^) indicates a negative aerodynamic damping which would cause a free-to-roll 

wing to be unstable. 

The pulse transfer-function analysis can be used to determine the rolling moment 

coefficient transfer-function and therefore the stability of the wing for a wide range of 

reduced frequency k (based on one-half of the wing root chord). A comparison of the 

pulse transfer-function analysis results, shown in Fig. 4, indicates the effects of angle 

of attack. At the lower angle of attack of a = 10° shown in Fig. 4(a), the lm(C, ) is 

*<t> 

negative for all values of reduced frequency which is indicative of stability in roll for 
small perturbations. In contrast, at the higher a = 30° angle of attack shown in Fig. 
4(b), the imaginary part is positive for k <. 0.5 which is indicative of instability in 
roll. Also, for a = 30°, the Re(C^) is negative in this range of reduced frequency 



which corresponds to a positive aerodynamic stiffness. The roll response will therefore 
oscillate (with increasing amplitude) rather than give rise to a static instability known 
as wing drop. 

The accuracy of the pulse analysis is verified by harmonic analyses performed at 
five values of reduced frequency: k = 0.0, 0.25, 0.50, 0.75 and 1.0. In these analyses, 
the wing was oscillated harmonically in roll with an amplitude of 1° for three cycles of 
motion with the rolling moment coefficient determined from the last cycle. The results 
of the harmonic analyses are compared with the rolling moment coefficient transfer- 
functions in Fig. 4. The good agreement between pulse and harmonic analyses therefore 
validates the accuracy of the pulse results for the cases considered. 

Forced Harmonic Results 

Because the pulse transfer-function analysis is limited to small perturbations, the 
large perturbation aerodynamic characteristics of the delta wing were investigated with 
a forced harmonic analysis. A reduced frequency of k = 0.25 was chosen for this 
analysis. This value lies at the midpoint of the range of reduced frequency identified by 

the pulse analysis as being an unstable condition for the free-to-roll wing at a= 30°. 
Three amplitudes of motion, <J) 0 = 5°, 15° and 35°, were considered at both a- 10° and 

a= 30°. The nondimensional timestep used for all cases was 0.00262. A comparison of 
rolling moment coefficient versus roll angle for each of these cases is shown in Fig. 5 to 
illustrate the effects of both roll amplitude and angle of attack. For the a= 10° angle of 
attack cases shown in Fig. 5(a), the results show a counter-clockwise-oriented loop for 
each roll amplitude which would produce a convergent (stable) response if the wing 
were free to roll. This prediction of a stable response at the smallest roll amplitude is 
consistent with the pulse transfer-function results of Fig. 4(a). Also as the roll 
amplitude is increased from <t> 0 = 5° to 15°, the rolling moment coefficient increases 

linearly. (Note the change in scaling of the vertical axis.) However, as the roll 
amplitude is further increased to <|) 0 = 35°, some nonlinear aerodynamic 


characteristics are exhibited in the "pinching" of the loop at the extreme roll angles 

although the free-to-roll response is still predicted to be stable. For the a= 30° angle 
of attack cases shown in Fig. 5(b), the results show clockwise-oriented loops for the <|> 0 

= 5° and 15° roll amplitudes which would produce a divergent (unstable) response if 
the wing were free to roll. This prediction of an unstable free-to-roll response at the 
smaller roll amplitudes is consistent with the pulse transfer-function results of Fig. 
4(b). At <J> 0 = 35°, counter-clockwise-oriented loops have formed at the extreme roll 

angles which consequently, would have a stabilizing effect on the free-to-roll response. 
The formation of these stabilizing loops was not, of course, predicted by the pulse 
analysis. In contrast with the oc= 10° angle of attack case, at the higher angle of attack of 
30°, the nonlinear aerodynamic effects at the larger roll amplitudes result in a change 
in the stability characteristics of the wing. 

Free-to-Roil Results 

The free-to-roll results were obtained for the flow conditions and structural and 
inertial parameter values listed in Table 1. The structural and inertial properties used 
in these calculations are loosely based on the characterisitcs of the models used in the 
experimental study of wing rock in Ref. 9. The initial angular velocity imposed on the 
wing was <(>'= 0.003 and the nondimensional timestep used was 0.004. The resulting 
roll angle response, shown in Fig. 6(a) for the a= 10° case, indicates that after the 
initial perturbation, the oscillatory response converges to its initial steady-state value. 
This stable free-to-roll response is consistent with the pulse and forced harmonic 
results of Figs. 4(a) and 5(a). The roll angle response, shown Fig. 6(b), for the a= 
30° case, indicates that initially the oscillatory response diverges for small values of 
roll amplitude which is consistent with the small amplitude pulse and harmonic results 
of Figs. 4(b) and 5(b). As the roll angle increases to around 35°, the rate of divergence 
decreases due to the stabilizing aerodynamics (counter-clockwise loops in the rolling 
moment coefficient at the extreme roll angles) shown in Fig. 5(b). Finally, the 



response reaches a maximum amplitude of motion at approximately 38° corresponding 
to a limit cycle. The reduced frequency of the limit cycle is k = 0.103. These results 
are similar in nature to those obtained by Arena and Nelson 1 1 in a low-speed 
experimental investigation of wing rock. The wing- rock time history from Ref. 11, 
shown in Fig. 7, was obtained for an 80° swept delta wing at 30° angle of attack. The 
reduced frequency of the wing rock in this case was approximately 0.125. Although the 
case considered in the present study is different from that of Ref. 1 1 (the data from Refs. 
8-12 are all for low speed flows whereas the conical Euler code is limited to supersonic 
freestream applications), the similarity between the two sets of results in Figs. 6(b) 
and 7 is noteworthy and gives credibility to the present calculations. 

Active. Roi l Suppressio n R esults 

An active rate-feedback control law was then implemented in an attempt to suppress 
the wing-rock motion. To determine an appropriate value for the gain, a stability 
analysis was derived by again recalling the nondimensional rolling equation of motion 
given by 


<> " = C i C / ( 16) 

Assuming that the rolling moment coefficient can be written as the superposition of the 
rolling moment coefficient transfer functions for <> and 8, 


C / + C / 6 5 0 9 ) 

Substituting the control law from Eq. (15) into Eq. (19) gives 

= + C/ # K V $' (20) 


Then for simple harmonic motion 



C, = {[Re(C /# ) - kM„K v Im(C /t )] 

+i [Im(C^ ) + kM„K v Re(C^ )]} <|> ( 2 1 ) 

On the right-hand side of Eq. (21), as before, the first term represents an aerodynamic 
stiffness, and the second term represents an aerodynamic damping. Therefore, 
stabilizing the motion of the wing requires that 

* m (^Atp ) + ^ M«> K v Re(C^ )<0 ( 22 ) 

Solving for the gain yields 


-1 ,m (C/ ) 

> 

v kM„ Re(C, s ) < 23 > 

A pulse analysis was performed to determine the rolling moment coefficient transfer 
function for 5 at a = 30° angle of attack (pulse amplitude of 5 0 = 1°). The transfer 

function from this analysis, shown in Fig. 8, indicates that for values of reduced 
frequency less than 0.5, the real part of is negative. Considering Eq. (23) and 

recalling from Fig. 4(b) that the values of lm(C/ ) for k<0.5 are positive, the value of 

$ 

K v must, therefore, be positive to suppress wing rock. The value of K v actually needs to 
be greater than that determined by evaluating the right-hand side of Eq. (23) in order to 
stabilize the wing since the above analysis assumes simple harmonic motion. For the 
flow conditions considered here, for example, the value for the gain that produces a 
neutrally stable (or simple harmonic) response is K v = 0.35. The free-to-roll 
analysis with active rate-feedback control was performed over a range of control gains: 
K v = 0.25, 0.40, and 0.50. The time-histories of the wing motion are shown in Fig. 9. 
The time-history for a control gain of Ky = 0.5 indicates a damped response as expected. 
Similarly, the response for control gain of K v = 0.4 is also damped although at a smaller 
rate than the K v = 0.5 case. The response of the wing for a control gain of K v = 0.25 
shown in Fig. 9 indicates that the response is no longer damped. However, a comparison 



with the results of Fig. 6(b) shows that the active rate-feedback control has a 
stabilizing effect on the free response at these conditions. 


Concluding Remarks 

A conical Euler code was developed to study unsteady vortex-dominated flows about 
rolling highly-swept delta wings, undergoing either forced or free-to-roll motions 
including active roll suppression. The flow solver of the code involved a multistage 
Runge-Kutta time-stepping scheme which uses a finite-volume spatial discretization of 
the Euler equations on an unstructured grid of triangles. The code allows for the 
additional analysis of the free-to-roll case, by including the rigid-body equation of 
motion for its simultaneous time integration with the governing flow equations. Results 
were presented for a 75° swept sharp-leading-edge delta wing at a freestream Mach 
number of 1 .2 and at a = 1 0° and 30° angle of attack. At the lower angle of attack of 
a = 10°, for example, a forced harmonic motion produced a rolling moment coefficient 
versus roll angle response that is counter-clockwise in orientation. This counter- 
clockwise response is indicative of stability in roll which was verified in a free-to-roll 
calculation. In contrast, at the higher angle of attack of a = 30°, the forced harmonic 
response was clockwise in orientation which indicates instability in roll. A free-to-roll 
calculation for this case produced an initially divergent response, but as the amplitude of 
motion grew with time the response transitioned to a limit cycle phenomenon known as 
wing rock. The wing rocking motion was subsequently suppressed, however, through the 
use of an active rate-feedback control law and antisymmetrically-deflected leading-edge 
flaps. 


A c kno w ledge m ent 

The authors would like to acknowledge Ken Morgan and Jaime Peraire of the Imperial 
College of Science, Technology and Medicine, London, England, for providing the 


20 



advancing front method grid generation program that was used to generate the grid in the 
present study. 


References 

1 High Angle of Attack Aerodynamics, AGARD-LS-121, December 1982; Lecture 
Series presented March 10-11, 1982 at NASA Langley Research Center, USA, March 
15-19, 1982 at von Karman Institute, Rhode-Saint-Genese, Belgium, and March 22- 
23, 1982 at DFVLR, Gottingen, Germany. 

2 Vortex Flow Aerodynamics, NASA CP-2416, July 1986; Proceedings of a 
conference sponsored by the NASA Langley Research Center, Hampton, Virginia, and the 
Air Force Wright Aeronautical Laboratories, Flight Dynamics Laboratory, Wright- 
Patterson Air Force Base, Ohio; held at NASA Langley Research Center, October 8-10, 
1985. 

Hamilton, W. T.: Manoeuvre Limitations of Combat Aircraft, AGARD-AR-155A, 
August 1979. 

^Hummel, D.: Documentation of Separated Flows for Computational Fluid Dynamics 
Validation, AGARD-CP-437, Vol. 2, December 1988, pp. p15-1 to p15-24. 

5 Kjelgaard, S. O.; and Sellers, W. L.: Detailed Flowfield Measurements Over a 75° 
Swept Delta Wing for Code Validation, AGARD-CP-437, Vol. 2, December 1988, pp. 
p10-1 to p10-14. 

^Squire, L. C.: Leading-Edge Separation and Cross Flow Shocks on Delta Wings, AIAA 
Journal . Vol. 23, March 1985, pp. 321-325. 

^Miller, D. S.; and Wood, R. M.: Lee-Side Flow Over Delta Wings at Supersonic 
Speeds, NASA TP-2430, June 1985. 

8 Nguyen, L. T.; Yip, L.; and Chambers, J. R.: Self-Induced Wing Rock of Slender 
Delta Wings, AIAA Paper No. 81-1883, August 1981. 

9 Levin, D.; and Katz, J.: Dynamic Load Measurements with Delta Wings Undergoing 
Self-Induced Roll Oscillations, Journal of Aircraft . Vol. 21, January 1984, pp. 30-36. 


21 



10 Jun, Y. W.; and Nelson, R. C.: Leading-Edge Vortex Dynamics on a Slender 

Oscillating Wing, Journal of Aircraft. Vol. 25, September 1988, pp. 815-819. 

1 1 Arena, A. S.; and Nelson, R. C.: The Effect of Asymmetric Vortex Wake 

Characteristics on a Slender Delta Wing Undergoing Wing Rock Motion, AIAA Paper No. 
89-3348, August 1989. 

12 Ng, T. T.; Malcolm, G. N.; and Lewis, L. C.: Flow Visualization Study of Delta Wings 
in Wing-Rock Motion, AIAA Paper No. 89-2187, August 1989. 

^Hoeijmakers, H. W. M.: Computational Vortex Flow Aerodynamics, AGARD-CP- 
342, Aerodynamics of Vortical Type Flows in Three Dimensions, July 1983, pp. 18-1 
to 18-35. 

"* 4 Newsome, R. W.; and Kandil, O. A.: Vortical Flow Aerodynamics - Physical Aspects 
and Numerical Simulation, AIAA Paper No. 87-0205, January 1987. 

15 Hsu C. H.; and Lan, C. E.: Theory of Wing Rock, Journal of Aircraft. Vol. 22, 
October 1985, pp. 920-924. 

16 Konstadinopoulos, P.; Mook, D. T.; and Nayfeh, A. H.: Subsonic Wing Rock of 
Slender Delta Wings, Journal of Aircraft . Vol. 22, March 1985, pp. 223-228. 

17 Elzebda, J. M.; Nayfeh, A. H.; and Mook, D. T.: Development of an Analytical Model 
of Wing Rock for Slender Delta Wings, Journal of Aircraft . Vol. 26, August 1989, pp. 
737-743. 

^Nayfeh, A. H.; Elzebda, J. M.; and Mook, D. T.: Analystical Study of the Subsonic 
Wing-Rock Phenomenon for Slender Delta Wings, Journal of Aircraft . Vol. 26, 
September 1989, pp. 805-809. 

19 Newsome, R. W.; and Thomas, J. L.: Computation of Leading-Edge Vortex Flows, 
Proceedings of the Conference on Vortical Flow Aerodynamics, NASA Langley Research 
Center, Hampton, Virginia, October 8-10, 1985. 

29 Kandil, O. A.; and Chuang, A.: Influence of Numerical Dissipation in Computing 
Supersonic Vortex-Dominated Flows, AIAA Paper No. 86-1073, May 1986. 


22 



Chakravarthy, S. R.; and Ota, D. K.: Numerical Issues in Computing Inviscid 
Supersonic Flow Over Conical Delta Wings, AIAA Paper No. 86-0440, January 1986. 

Murman, E. M., Powell, K. G.; Miller, D. S.; and Wood, R. M.: Comparisons of 
Computations and Experimental Data for Leading Edge Vortices - Effects of Yaw and 
Vortex Flaps, AIAA Paper No. 86-0439, January 1986. 

23 Murman, E. M.; and Rizzi, A.: Applications of Euler Equations to Sharp Edge Delta 
Wings With Leading Edge Vortices, AGARD-CP-412, November 1986. 

24 McMillin, S. N.; Thomas, J. L.; and Murman, E. M.: Euler and Navier-Stokes 
Solutions for the Leeside Flow Over Delta Wings at Supersonic Speeds, AIAA Paper No. 
87-2270, July 1987. 

25 Fujii, K.: A Method to Increase the Accuracy of Vortical Flow Simulations, AIAA 
Paper No. 88-2562, July 1988. 

26Powell, K. G.; Murman, E. M.; Perez, E. S.; and Baron, J. R. : Total Pressure Loss 
in Vortical Solutions of the Conical Euler Equations, AIAA Journal . Vol. 25, 1987. 

27 Kandil, O. A.; and Chuang, A. H,: Computation of Steady and Unsteady Vortex 
Dominated Flows, AIAA Paper No. 87-1462, June 1987. 

28 Kandil, O. A.; and Chuang, A. H.: Unsteady Navier-Stokes Computations Past 
Oscillating Delta Wing at High Incidence, AIAA Paper No. 89-0081, January 1989. 

29Batina, J. T.: Vortex-Dominated Conical-Flow Computations Using Unstructured 

Adaptively-Refined Meshes, AIAA Paper No. 89-1816, June 1989. 

3 °Lee, E. M.; and Batina, J. T.: Conical Euler Solution for a Highly-Swept Delta 

Wing Undergoing Wing-Rock Motion, NASA TM 102609, March 1990. 

3 ^ Karagonnis, T.; Maxworthy, T.; and Spedding, G. R.: Generation and Control of 
Separated Vortices Over a Delta Wing by Means of Leading Edge Flaps, AIAA Paper No. 
89-0997, March 1989. 

32 Ng, T. T.: On Leading Edge Vortex and Its Control, AIAA Paper No. 89-3346, 
August 1989. 


23 



33wood, N. J.; Roberts, L.; and Celik, Z.: The Control of Asymmetric Vortical Flows 
Over Delta Wings at High Angles of Attack, AIAA Paper No. 89-3347, August 1989. 

34 Ellis, D. G.; and Stollery, J. L.: The Behavior and Performance of Leading-Edge 
Vortex Flaps, ICAS-88-4.5.2, 1988. 

35Batina, J. T.: Unsteady Euler Airfoil Solutions Using Unstructured Dynamic 

Meshes, AIAA Paper No. 89-0115, January 1989. 

36seidel, D. A.; Bennett, R. M.; and Whitlow, W., Jr.: An Exploratory Study of Finite 
Difference Grids for Transonic Unsteady Aerodynamics, AIAA Paper No. 83-0503, 
January 1983. 

37seidel, D. A.; Bennett, R. M.; and Ricketts, R. H.: Some Recent Applications of 
XTRAN3S, AIAA Paper No. 83-1811, July 1983. 

38Morgan, K.; and Peraire, J.: Finite Element Methods for Compressible Flow, Von 
Karmen Institute for Fluid Dynamics Lecture Series 1987-04, Computational Fluid 
Dynamics, March 2-6, 1987. 



Table 1 Summary of structural parameter values and flow conditions 
for the free-to-roll calculation. 



Fig. 1 Partial view of unstructured grid about a 75° swept delta wing. 
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(a) 8 = 10°. 



Fig. 2 Partial view of deforming mesh about deflected leading-edge flap. 


26 







(a) a = 10°. 



(b) a = 30°. 

Fig. 4 Angle of attack effects on the rolling moment coefficient transfer-function versus 
reduced frequency for a 75° swept delta wing at = 1.2 (pulse amplitude 

<D 0 = n- 


28 




- 0 . 02 |_ 

I I 

— 10 -5 

o.o6r 


j i i 

0 5 10 

<t> 0 = 15° 


C t 0.00 



-0.06L 

I I L 

-30-15 0 
0. 1 4f 


i i 

15 30 

♦o = 35° 


I I L 

-10 -5 0 


_L 

5 


10 




-0.14L 

I L 


-70-35 


_L 

0 


_l I 

35 70 


i 1 i i i 

-70-35 0 35 70 


<l> (deg.) 


4> (deg.) 


(a) a = 10°. 


( b ) a = 30°. 


Fig. 5 Angle of attack effects on the rolling moment coefficient versus instantaneous roll 
angle for a 75° swept delta wing at = 1 .2 and k = 0.25. 
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Fig. 7 Wing-rock time history for an 80° swept delta wing at 30° angle of attack 
(Ref. 11, reprinted with permission from Professor Robert C. Nelson 
Notre Dame University). 



k 


Fig. 8 Rolling moment coefficient transfer-function for the flap versus reduced 

frequency for a 75° swept delta wing at =1.2 (pulse amplitude 5 0 = 1°). 
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Fig. 9 Free-to-roll response with rate-feedback control for a 75° swept delta wing at 
= 1.2 and a = 30°. 
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